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ABSTRACT 


The achievable gain in the radiation resistance and directivity of a 
low frequency underwater transducer due to the presence of an array 
of sympathetic resonators was analyzed. The resonators were all taken 
to be air bubbles, and both the resonators and transducer were taken 
to be compact (ka<<l). The resonators were taken to be equally 
spaced around a circle of radius R, with the transducer located on the 
axis. The gain was calculated for various numbers of resonators as a 
function of Kayesonator, Katransducer and kR, for the transducer in the 
plane of the resonators and out of the plane a distance of one-quarter 
wavelength. For the transducer in the plane, a gain in radiation 
resistance of approximately two is possible with six or more 
resonators. For the transducer out of the plane, it is shown that a 
significant gain in directivity can be achieved at the expense of a slight 


decrease in the gain in the radiation resistance. 
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THESIS DISCLAIMER 


The views expressed in this thesis are those of the author and do 
not reflect the official policy of cr position of the Department of 
Defense or the U.S. Government. 

The reader is cautioned that computer programs developed in 
this research may not have been excerised for all cases of interest. 
While every effort has been made, within the time available, to ensure 
that the programs are free of computational and logic errors, they 
cannot be considered validated. Any application of these programs 


without additional verification is at the risk of the user. 
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LIST OF SYMBOLS 


ai Radius of element n. 
¢ Speed of sound in water. 
Gr The ratio of the acoustic radiation resistance of a 


transducer with sympathetic resonators present to the 
acoustic radiation resistance of an identical transducer 
with no resonators. 


k Wave number, wi. 
ka, The resonance value of ka for a single air filled bubble. 
lan Separation distance between element n and m. 
Ph Incident free field acoustic pressure at element n 
due to external field only. 
Po Surface pressure at element n. 
Pre Scattered pressure at the surface of resonator n, in the 
presence of other rescnators. 
Poingh. Scattered pressure at the surface of a single resonator. i.e. 
in the absence of other resonators. 
G Mass injection rate of a spherical sound source. 
r Dis anc, from scatterer center to a field point. 
Z Specific acoustic impedance, pressure divided by surface 
normal velocity. 
A Single scatter amplification coefficient of element n. 
Multiple scatter amplification coefficient of element n. 
Cam Ratio of acoustic impedances defined as: 
— Ln 
caer 
ix 


Fy 


a 
Zn 

at 
Zan 

ar 
Zam 


at 
total 


FAc° 





Ratio of acoustic impedances, defined as: 


pace 
a 


Scattered acoustic pressure from a single scatterer, 
normalized to the incident acoustic field. 


Scattered acoustic pressure from a scatterer,interacting 
with other scatterers. normalized to the incident acoustic 
field 


Electrical current. 

Ambient pressure of environment. 
Surface area of element n. 
Transduction coefficient. 

Volume velocity of element n. 


Volume velocity of element n, normalized to transducer's 
volume velocity. 


Volume of source. 
Acoustic impedance, pressure divided bvolume velocity. 
Acoustic impedance of element n. 


Acoustic radiation impedance seen by element n with all 
other elements blocked. 


Acoustic transfer impedance from element m to element 
n. 


The total acoustic radiation impedance seen by transducer 
in the presence of sympathetic resonators. 


Transducer open circuit (1=0) acoustical impedance. 


Oni 


blo 


Kronecker delta function, where: 


San _| ee | 
| 0; m#n | 
Densitv of water. 


Ratio of the magnitudes of the scattered acoustic pressure 
at the surface of a resonant scatterer to the incident 
acoustic pressure in the presence of other scatterers. 


Ratio of the magnitudes of the scattered acoustic pressure 


at the surface of a single resonant scattererto the incident 
acoustic pressure. 


Ratio of specific heats (air). 

Angular frequency of plane wave. 

Angular frequency of scatterer resonance. 
Defined as: 


leojyen? 1] 


t= 
7 ka 


A Nx] column matrix, n = 1...N. 


A NxN square matrix, n,m = 1...N. 
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I. INTRODUCTION 


A. BACKGROUND AND HISTORY 

I]. Tolstoy [Ref.1],[Ref 2] and A. Tolstoy and I. Tolstoy [Ref. 3] nave 
recently analyzed the multiple scattering of sound from a collection of 
compact (ka<<1) monopole scatterers (e.g. gas filled bubbles, gas filled 
balloons and thin shells in water) insonified by a simple harmonic 
plane wave at frequencies close to the resonance frequency of an 
individual scatterer. The response of these systems was quantified by 
computing the ratio of the scattered pressure on the surface of a 
scatterer in the presence of the other elements to that of the incident 
plane wave. This ratio was defined as the amplification factor, p. They 
investigated the value of this amplification factor for various simple 
geometrical arrangements of monopole resonant scatterers, and 
showed that amplification factors of the order of 107 were to be 
expected for air bubbles in water at one atmosphere ambient pressure. 
(The amplification factor of a single bubble is about 70.) These 
systems were said to be "quasiresonant”. 

I. Tolstoy [Ref. 4] has also recently analyzed the strength and 
directionality of the resonant multiple scattering of sound from a 
linear array of compact monopole resonant scatterers irradiated by a 
plane wave. He reports that "quasiresonance” effects exhibit strong 
directionality in this system, with amplification factors of up to 700. 


The potential application of resonant multiple scattering, or 


"quasiresonance”, in Tolstoy's words, to improve the performance of a 








low frequency active sonar transducer is the subject of the research 


described in this thesis. 


B. OBJECTIVES 

The objectives of this research are multiple. The first is to 
reproduce I. Tolstoy's results using network analysis, a tool familiar to 
the transducer and array designer. It will then be shown that the two 
methods of analysis, multiple scattering and network analysis, give 
identical results. The existence of quasiresonance will then applied to 
a circular array of scatterers, where the incident plane wave is 
replaced by an active transducer located on the axis of the circle and 
radiating spherical waves. The acoustic advantage gained in this 
system will be quantified by the gain in the radiation resistance seen 
by the transducer in the presence of these scatterers. An increase in 
radiation resistance enables the transfer of more power to the acoustic 
field for a given surface velocity of the transducer. (High power low 
frequency active sonar transducers tend to be resonant devices that 
are displacement limited.) The cases of the transducer in the plane of 
the circle and displaced by one-quarter wavelength from the plane will 
be analyzed for both their increase in radiation resistance and gain in 


directivity over that of a single transducer. 


Il. THEORY 


A. QUASIRESONANCE 

For simple configurations of compact resonant scatterers 
insonified by a simple harmonic plane wave, Tolstoy predicted the 
existence of so-called "quasiresonance" using multiple scattering 
theory. (Quasiresonance is Tolstoy's term for the enhancement of the 
scattered pressure from a resonant scatterer as a result of multiple 
scattering.) Tolstoy's analysis will be reviewed and applied to the 
simplest case in section Al. In section A2 the same problem is treated 
using network analysis. The existence of quasiresonance is again 
predicted using a coupled network model. It will be shown that this 


second method produces results identical to Tolstoy's. 


1. TOLSTOY'S ANALYSIS OF MULTIPLE SCATTER 

Tolstoy's analysis of the multiple scattering of a plane wave 
incident upon a collection of compact resonant scatterers proceeds 
using linear theory of acoustics. The following development parallels 
that of Reference [1] and Reference [3], with the difference being that 
here all time-harmonic functions are taken to be of the form et. This 
will result in minor differences in the appearance of otherwise 
equivalent equations. 

Consider first the case of a single resonant scatterer insonified 
by an incident plane wave. From Lighthill [Ref. 5], the sound field of a 


spherical source is given by Equation 2-1. 


plr.t (4nry q(t-1/c) (2-1) 
where: 
r = Distance from scatterer center to the field point. 
q = Mass injection rate of a spherical sound source. 
c = Speed of sound in water. 


The dot over q indicates the derivative with respect to time. Treating 
the single scatterer as an air fille. pubble, for an incident plane wave 


of free field pressure pf, q takes the form of Equation 2-2. 





P >,,{¥P . 
4 = pray (le : i} ~(1-jka) 2 | eicat 
pe? ws (2-2) 
where: 

p! = Incident free field acoustic pressure. 
k = Wave number, w<. 
V = Volume of scatterer. 
a = Scatterer radius. 
y = Ratio of specific heats (air). 
c = Speed of sound in water. . 
p = Density of water. 
Po = Ambient pressure of environment. 
w = Angular frequency of plane wave. 
Wo = Resonance frequency of scatterer. 


Inserting Equation 2-2 into Equation 2-1, dropping the time 
dependence, and solving for the ratio of the scattered acoustic 
pressure at some radius r from the scatterer, psingie(r), to that of the 


incident pressure, p' , one obtains Equation 2-3. 


Pringlel) Pe 





e-iks 


poo jkr (el (2-3) ; 





where: 


For situations where the surface tension and corrections due to energy 
losses in the interior of the scatterer (e.g. bubbles, balloons) are 
negligible, the ratio of the resonance angular frequency of the 
scatterer to the angular frequency of the incident plane wave is given 


by Equation 2-4. 


Oo _ / 3yPo 


o V pc? (ka? (2-4) 





Letting r equal a, the scatterer radius in Equation 2-3, results in 
Equation 2-5, the ratio of the scattered acoustic pressure at the single 


scatterer'’s surface to that of the incident free field pressure. 


P ingle = i 


po ka 


1 


1-j6 


evjka 








(2-5) 


Tolstoy introduces the quantities F, and A by rewriting Equation 


2-3 as Equation 2-6. 


p, Painetelt) _ 4 evikr 
pf jkr (2-6) 





where: 
A= 
1-j¢ 


Fg is the normalized scattered acoustic pressure for a single scatterer 
insonified by an incident plane wave and A is the single-scatterer 
amplitude coefficient. 

Now consider a system of N identical interacting scatterers 
insonified by a plane wave. Let p;‘(r) represent the scattered acoustic 
pressure. An equation for p&(r)/pi,, equivalent to that of Equation 2-6 
for a single scatterer, can be written as Equation 2-7 for the multiple 


scatterer situation. 


pf j (2-7) 


Fp, is the normalized scattered acoustic pressure from scatterer n, 


insonified by both the incident plane wave and the fields of the other 
scatterers, Pa is the free field pressure of the incident plane wave at 
scatterer n, and By is the multiple scatterer amplitude coefficient for 
scatterer n. 

For compact scatterers, the overall response of a scatterer to 
multiple incident pressure fields can be evaluated by summing up the 
response to each individual incident pressure field [Ref. 3]. For a 
system of N scatterers, this results in Equation 2-8 for the scattered 
field at the surface of the nth resonator. 

B, f=" = Ag (1+ Y) Bm S27 civ | ese 


a jka, kl m jka, 
J <a Oe J (2-8) 








where: V_., is the phase difference in the incident plane wave between 


resonator m and n, and m=1...N. The subscript notation for A has 
been added to allow for a system of different sized scatterers. 

Tolstoy defines the amplification factor p (different than the 
scatter amplitude coefficient mentioned previously) to be the 
magnitude of the ratio of the scattered acoustic pressure on the 
surface of a scatterer, p;, to the incident free field acoustic pressure 
of the incident plane wave, ph. With this definition, Equation 2-9 


follows from Equation 2-7. 


~ Bal 
ka, (2-9) 





Le 
pf 


Tolstoy denotes the amplification factor for a singic scatterer as Uo. 


Equation 2-10, then. follows from Equation 2-6. 


pf Kan (2-10) 


From Equations 2-9 and 2-10, pt/vo is expressed in terms of 


Jor” fae |= Boda 
Ho, singlep (2-1 1) 


a. APPLICATION OF TOLSTOY'S METHOD 


Bn/An by Equation 2-11. 


As an example application, Tolstoy considered a system of two 
identical resonators insonified by a plane wave traveling along the axis 


joining their centers, as shown in Figure 2-1. 


ee 





: Resonator of 
Incoming Plane Wave radius =a 


/— l —>| 
Figure 2-1. A Two Resonator System Insonified by 


a Plane Wave Along the Axis of The Resonators. 


Application of Equation 2-8 to the situation shown in Figure 2-1 


results in the set of Equations 2-12. 


B e-jka aaewe sok evjka Bs e7 ik! e-ikl 


' jka ~~" jka ja Kl 








By ei? aac 4 ae pe citi 
jka ka jka kl (2-12) 


For a system of identical scatterers, the values of A; and Ag are equal. 
Therefore, the subscript on A in Equation 2-12 has been omitted. 
Since all variables in Equation 2-12 are known except B; and Bg, one 


can solve for B)/A in terms of A and kl, resulting in Equation 2-13. 


1+ [A | eK 
J 


B,/A= 
Aci? 


1- id 








(2-13) 


Substituting the expressions for (and /wp for an air-filled bubble into 


the expression for A, B;/A is given in terms of the variables ka, kl, c, 


Po, and p by Equation 2-14. 





ka gi(ka-k1) 
1 - e-ik! ki 




















pc (kar) Ria ee 
By AS PK ae pc*{ka) 
Psingle ea 2 (2- 14) 
pc (ka ka ee 
| pc*{ka) 











For bubbles just under the surface of the ocean, we take: 


Po = 10132 Pa (1 atmosphere) 
p = 998 kgm 
c = 1498 ms:! 


As observed by Tolstoy [Ref. 3], plots of Equation 2-14 display 
spurious peaks. Following Tolstoy, then, the graphs of !B,/Al versus 
ka and kl, plotted in Figures 2-2a,b have been treated to limit peaks in 
bound intervals less than 2ka. (According to Tolstoy peaks narrower 
than 2ka, where 2a equals the resonator diameter, are not meaningful, 
since the pressure over a compact resonator is approximately 
constant.) The truncation method used involvec comparing three 
consecutive values of !B;/A! in the kl direction for each value of ka. 
(Three consecutive points in the kl-direction span approximately 
2ka.) If the value of 1B,/Al for the middle point was greater than 
either of its neighboring points (a peak) then this value was set equal 
to the mean of the two neighboring points. This procedure was 


repeated several times to get the desired smoothing. 


Fig. 2-2a,b 
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. ka 
- Three Dimensional (a), and Contour (b) Plots Showing 


1B,/Al for the Double Resonator System as a Function of ka and 
kl for Axial Incidence of the Plane Wave. 








Figures 2-2a,b show a maximum response of 1B,/AI! equal to 5.5 


at approximately ka=0.01393 and kl=0.60. As a result of multiple 
scattering, then, the scattered pressure amplitude of the first bubble 
is 5.5 times greater in the presence of the second bubble than in its 
absence. Additionally, the value of ka at the maximum response is 


shifted from that for a single resonant bubble of .01379. 


2. NETWORK ANALYSIS OF MULTIPLE SCATTER 
Fig. 2-3 displays a network model of N interacting scatterers 
insonified by an incident acoustic wave. (Note: pressure (p) and 


volume velocity (Un) are chosen as the mechanical variables, and the 
direction of positive velocity is clockwise.) The network equations for 


Fig. 2-3 are: 


N 
Pn = Ph + > ZamU m Pn = -Zn Un (2-15) 


m=1 


where for Fig 2-3: 


Pn= Surface pressure at element n. 


pis Free field pressure at element n due to 
an externally incident field only. 


Un= Volume velocity of element n. 
= Acoustic impedance of element n. 


n= Acoustic radiation impedance seen by 
element n with all other elements 
blocked. (Umzn= 0 me=1...N) 


= Acoustic transfer impedance from 
element m to element n. 


11 





Resonators Acoustic Field 





Figure 2-3. Network Model of N Interacting 
Resonators Insonified by a Plane wave. 
Neglecting the effects of surface tension, thermal relaxation, 
etc., the intrinsic acoustic impedance of a small (ka<<1) gas-filled 
bubble can be approximated by one of reactance only, given bv Kinsler, 


Frey, Coppens and Sanders [Ref. 6] as Equation 2-16. 


-3yP 
Piaj Sh 
Akan) (2-16) 
where: 
a, = Radius of element n. 


Speed of sound in water. 
Wave number, «yc. 
Py = Ambient pressure cf environment. 


2 





Sn Surface area of element n. 


Ratio of specific heats (air). 


The self acoustic radiation impedance seen by element n, Zan. 


can be approximated by Equation 2-17. 


Zit, = pc ((kayF+jkay] S; (2-17) 
where: 
P = density of water. 


Equation 2-17 is the acoustic radiation impedance presented to a 
compact source with no other active or scattering elements present. 
It is also the acoustic radiation impedance seen by a compact source in 
the presence of other compact Sources when these sources are 
blocked (U,,4, = 0 me=1...N). 
The acoustic transfer impedance from element m to element n 
is given by Equation 2-18. 
Zit = ZA Sagi kan-klon) 
Klam (2-18) 
In matrix form the network equations (Equation 2-15) can be 


written as Equation 2-19. 


tint | Zeal Um + ph ana es (2-19) 
where: 
eae n=1...N A 1xN column matrix, N elements. 
‘Xam nam =1..N A NxN square matrix, N? elements. 


Substituting py = -Z;,Un and rearranging results in Equation 2-20. 
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lpf!= [Z3+Z2Sam| {Um}  m,n=1...N (2-20) 


where 
Sam = Kronecker delta function, 


5 2 1; m=n | 


| 0; men [, 


a. APPLICATION OF NETWORK ANALYSIS 
Modeling the system of two scatterers shown in Figure 2-1 as a 


network, Equation 2-20 becomes Equation 2-21. 


M1tZ, Zi) 





U; | 
U2 














21 2% +Z} (2-21) 
Using the relation p§° = ZjnUn_ to eliminate Un, the above equations 
can be written in terms of the scattered acoustic pressure at each - 
individual resonator, pn, and the incident plane wave pressures as 
given by Equation 2-22. 
pit = (pi + Cio pF Dy 





px = (pf + Car pit ) DY (2-22) 
where 
Se en eee 
Zam Zan 


For a system of identical resonators, Cp, and Dy, are the same for all m 
and n; therefore the subscript notation is not required. Solving 


Equation 2-22 for the scattered acoustic pressure at each resonator in 





terms of the incident plane wave pressure results in Equation 2-23 


for the scattered pressure at resonator one. 


c Pi D+ppC 


Pi 
D2 -C (2-23) 


The scattered acoustic pressure at the first resonator in the 
absence of the second resonator, Psingle;, is obtained using Equation 2- 
23 by letting the transfer coefficient C equal zero, resulting in 
Equation 2-24. 


_Ph 


Psingicy D. (2-24) 


Dividing equation Equation 2-23 by 2-24 results in Equation 2- 
25 for the ratio B;/A as defined by Tolstoy. 





/A= = 


Ponple (cy 
wey? Pele (2-25) 


Inserting Equations 2-16 thru 2-18 into Equation 2-25, one 


obtains an equation for B,/A in terms of the variables ka, kl. c, Po, pi i 


po and P as given by Equation 2-26. 
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By,/A=—Pi_= pc%{ka) 
Psingle kaestke-t1) 2 (2-26) 
He diesen. Sayre | 
pc (eas ka - —~— ; 
pcka) I} , 








The incident plane wave pressures Pi and pe are of equal magnitude 


and differ only by the phase factor eKz, resulting in Equation 2-27. 


a esikl 


| 
= roe 


(2-27) 


Substituting this expression into Equation 2-26, one obtains Equation 




















2-28. 

ka oj(ka-k1) 

ne 3yP 

pe [taj ka - SO : 
pir pc%ka) ff | 
pee 2 (2-28) 
Pringle ka gj(ka-kl) - 
: kl 
F 3Po 
c{(ka}?+jlka - 
3 | : pc%ka} 











Equation 2-28 is identical to Equation 2-14, demonstrating the 


equivalence of Tolstoy's method and the network analysis method. 


3. A TRANSDUCER IN A MULTI-RESONATOR FIELD 
Network analysis will now be applied to the problem of a compact 
transducer in the presence of an array of interacting, compact 


resonant scatterers. The procedure is the same here as in section A2, 
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except that the incident plane wave is replaced by an outward 
traveling spherical wave whose source is a compact transducer. A 


network model of the transducer and N resonator system is shown in 


Figure 2-4. 


Acoustic Field 


Transducer 


Resonators 





Fig. 2-4. Network Model of a Single Transducer 
and an N Resonator Field. 


The network equations for Fig. 2-4 are; 
Po= >, ZanUm pp =-Z4 Un form=1..N — po=TI-Z* Up (2-29) 
m 


where 
T 
hin 


transduction coefficient. 
transducer open circuit (I=0) acoustic impedance. 


Ws 





and the subscript n=O denotes the transducer, n#0 a resonator. 

The transducer's open circuit acoustic impedance (Z*°) and 
transduction coefficient (T) are shown only for completeness and will 
not be considered further, as the quantity of interest will be the 
radiation impedance presented to the transducer. 

Using the network Equations 2-29, the matrix equation for the 


network model shown in Figure 2-4 is given by Equation 2-30. 














; T 
Po Z00 [Zao Uo 
atl gear peach ale aT te 
nme {Znot (Zam) [Heme (2-30) 
where 
Xno n=1...N denotes a 1xN column matrix, N elements. 


f \ 

\ f 

[Xnm!  mn,me=l1...N denotes a NxN square matrix, N2 elements. 
and the superscript T denotes transpose. Replacing the 1 x N surface 
pressure matrix, {pn}, with {-Z,Un} and then normalizing by Uo, the 


transducer's volume velocity, Equation 2-30 becomes Equation 2-31. 

















total | 00 [Ze ol" 1 
= >; nem=l..N 
10} (Zno} [ am +Z38nm| lure! (2-31) 
where 
total = Total acoustic radiation impedance seen by the 
transducer. 
UR! = Ratio of element m’'s volume velocity to that of the 


transducer (US! = U,/Uo)- 
Equation 2-31 can be written as two separate matrix equations, 


Equations 2-32a.b. 
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Ze a= Zoo HZ AU (2-32a) 


(0} = {Z| + [Zah, +Z85aml{ UES} (2-32b) 


Solving Equation 2-32b for (U%S!} by matrix inversion and substituting 
the result into Equation 2-32a, one obtains Equation 2-33, the total 


acoustic radiation impedance seen by the transducer. 
T “1 
Zit sa = Zi ~ (Zatg\” [ZH 4Z28 uml (Z%} (2-33) 


Equation 2-33 expresses this impedance in terms cf the radiation and 
transfer impedances of the individual elements. 

The tot*l acoustic power radiated by the transducer is given by 
Equation 2-34. 


I] = | Ud? Re(Ztora } (2-34) 


The gain in power delivered to the acoustic field due to the addition of 
the sympathetic resonators equals the ratio of total power delivered 
with the resonators present to the power delivered when they are 
absent, 1 = 1] Up? Re(Ziy ). For a displacement-limited source, typical of a 
sonar projector operated near its resonant frequency, this ratio 
reduces to the ratio of the radiation resistances for the two cases. This 


gain factor is expressed as Equation 2-35. 


~ RelZ¢, 


“total / 


Re(Zg } (2-35) 
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It. ANALYSIS 


A. ANALYSIS OF RESONATOR AND TRANSDUCER CONFIGURATIONS $ 
The theory developed in the previous chapter will now be applied 

to two specific geometric configurations of a transducer and N 

identical resonators. The resonators will be treated as simple gas filled 

bubbles. It is assumed throughout that the acoustic wavelength is 

much larger than the dimensions of the transducer and the 

resonators. For initial calculations, ambient pressure is taken to be one 

atmosphere, and the radius of the transducer is taken to be equal to 


that of the resonators (kao = ka,). 


1. PLANAR CONFIGURATION 
The geometry of the first arrangement to be analyzed is shown 


in Figure 3-1. It consists of N identical resonators equally spaced 





Figure 3-1. Resonators Equally Spaced around Transducer 
in a Planar Configuration of Radius Ino, N=4. 








around a circle of radius lho. This geometric arrangement will be 


referred to as the planar configuration. 
For the configuration shown in Figure 3-1, Equation 2-35 


becomes Equation 3-1. 


-1 
Rel Zoo 7 (Zeer [z2,, +Z38 nam! (Zao! n,m=1...4 


Cee 
: Re{ Zoo) (3-1) 





Using Equations 2-16, 17, and 18, and canceling the common surface 


areas S,, Equation 3-1 can be written as Equation 3-2 in terms of 


specific acoustic impedances: 


Gp = Rel ato t2tfol™ [atin +288am! tet!) ag 
Rel 2) (3-2) 
where: 
Zit) = pe ((kay-rjxa) 
Z2t, +285 am = pc ((ka}*+jka) - j ao ; when n=m 
Z8,, +225am = Zam = PC ((ka)2+jka) a <A_ e}(ka-klnm) when n#m 


KIlnm 
zy = pc (ka e+jka) _ka_ ej (ka-kl,,9) 
klno 


and Inm is the distance between element n and m, the subscript 0 
denoting the transducer. Using simple trigonometry, kKlhm for n,m #0 
can be written in terms of m, n, and klho as Equation 3-3 for the 


planar case. 


Klam = kyo 2 sin In-m} where nm = 1...N (3-3) 
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Using Equations 3-2 and 3-3, the gain in radiation resistance seen by 
the transducer, Gr, can be expressed in terms of ka and Klyo. 

The gain in radiation resistance for the N=4, planar 
configuration can be plotted in the form of surface and contour plots 
using Equation 3-1, as shown in Figures 3-2a and 2b. It should be 
pointed out that here and in all subsequent surface plots of Gr, no 
spurious peaks were observed, and no smoothing has been applied. 
Inspection of these graphs reveals two local maxima, one of magnitude 
4.34 at ka=0.975 ka,, kl,o=0.5 and another of magnitude 1.79 at 
ka=1.005 ka,, klho=2.7, where ka, equals 0.01379, the resonance 
value of ka for a single air filled bubble at one atmosphere ambient 
pressure. The local maximum at ka=0.975 kay,, kl,o=0.5 is an artifact 
of the range of kljg chosen. In fact, the gain in radiation 
resistance diverges as (klno) when klpo approaches zero. This case 
is uninteresting and will not be considered further. The interesting 
local maximum at ka=1.005 ka,, klho=2.7 represents a gain in the 
radiation resistance seen by the transducer by a factor of almost two. 
Note that kl = x corresponds to a separation distance of one-half 
wavelength. A maximum gain in radiation resistance is therefore 
obtained when the resonators are located a distance just under one- 
half wavelength from the transducer. 

The gain in radiation resistance for the planar configuration 


where N is set equal to eight and fifteen is shown in Figures 3-2c thru 


3-2f. These plots display similar characteristics to the N = 4 case. All 





















show a local maximum of magnitude approximately two at ka = ka, = 
0.01379 and klho =2 

The local maximum values for Gg at Klho =n, ka = ka, = 
0.01379 for various values of N ranging from 2 to 25 and their 
associated parameters are tabulated in Table 3-1. Two major points 
can be observed from Table 3-1. First, the resonant bubbles’ volume 
velocity leads the transducer's by a phase shift of approximately 90° 
(equivalent to -270°). This phase shift is created in two parts. Since 
the separation between transducer and resonator is approximately one 
half a wavelength, this introduces a_ shift of -180° in phase due to 
simple wave propagation. The remaining -90° is introduced at the 
surface of a compact resonating bubble between the incident and 
scattered wave as predicted by Equation 2-24. Second, note that little 


is gained as one increases the number of resonators beyond about four. 





Table 3-1. Values of Gp for Values of N Ranging from N=2 
to 25 andthe Parameters Associated with the Local 
Maximum at kl, =n, ka = ka,;=0.01379 for the Planar 

Configuration. 
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Figure 3-2a,b. Three Dimensional and Contour Plots of Gras a 
Function of ka and kino for N= 4, Planar Configuration, 1 atm. 
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3 Figure 3-2c,d. Three Dimensional and Contour Plots of Gras a 
Function of ka and klpo for N=8, Planar Configuration, 1 atm. 
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Figure 3-2e,f. Three Dimensional and Contour Plots of Gras a 
Function of ka and klpo for N=15, Planar Configuration, 1 atm. 
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a. PLANAR CONFIGURATION BEAM PATTERN 

Polar plots of the beam patterns can be developed for the planar 
configuration using the axis system defined in Figure 3-3, The 
origin is located at the transducer, with the resonators located in the 
x,y plane and 6 being the angle measured from the positive y axis in 
the x,y plane. 

Figure 3-4 clearly demonstrates the acoustic advantage of 
the planar configuration over that of a lone transducer. The zero dB 
reference is taken as the field of an identical transducer with equal 


volume velocity in absence of the resonators. 





x 
y.z plane 7.9 =45 ‘plane 


Figure 3-3. Axis System for Beam Patterns Plots, Planar 
Configuration. Note the definitions of y.z and z, 0=45° Plane 
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Figure 3-4. Beam Pattern for the Planar 


Configuration in the y,z Plane and the z,@=45° 
Plane, Po= latm. Scale in Decibels (dB). 


2. CONICAL CONFIGURATION 

The second geometric arrangement analyzed was chosen to 
take advantage of the 90° phase shift between the volume velocity of 
the resonators and the transducer for maximum radiated power. If the 
transducer is offset from the plane of the resonators by one quarter 
wavelength along the positive z axis, the radiation from the resonators 
will constructively interfere with that of the transducer along the 
positive z axis. This configuration will increase the directivity of the 


system. This arrangement will be referred to as the conical 


configuration and is shown in Fig. 3-5. 











Transducer £22:) — — —- — — 


Resonator n 


Figure 3-5. Resonators Equally Spaced Around Transducer 
in a Conical Configuration at Radius Rand N=4. The 
Transducer is Offset From the Plane of the Resonators by a 


Distance of 4/4. 


The value of Gp can again be expressed in terms of ka and klpo 
using Equation 3-1, where Equation 3-3 is used to express kKlnm in 
terms of n, m, and kl,o for the conical configuration. 


Klam = a] tklno? {8} 2 sin{ * |n-m] where n,m = 1...N (3-3) 


The gain in radiation resistance for the conical configuration for 
the cases N=4, 8 and 15 are plotted in the form of surface and 
contour plots, as shown in Figures 3-6a through 6f. Inspection of these 
plots reveals the existence of the same local maximum as the planar 
configuration of magnitude about 2 at klho=n, and ka = .01379. Note 
that, unlike the planar case, the kl axis of Figures 3-6a through f 
begins at klno=1.75 vice 0.5. This is due to the displacement of the 
transducer in the positive z direction. This displacement requires klno 


be at least 1/2. 
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The local maximum values for Gp at kly»o = 7, ka = .01379 for 
various values of N ranging from 2 to 25 and their associated 
parameters are tabulated in Table 3-2 for the conical configuration. In 
each case the value of Gp for the conical configuration is slightly less 
than for the planar configuration. Also, the value of klmax is slightly 
greater for the the conical configuration, but is still approximately 
equal to xn. (A value of x for klno corresponds to a value of 30 degrees 


for the angle a in Figure 3-5.) 





Table 3-2. Values of Gr for Values of N Ranging from N=2 
to 25 and the Parameters Associated with the Local Maximum 


at kl=n, ka = kl, = 0.01379 for the Conical Configuration. 
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: Figure 3-6a,b. Three Dimensional and Contour Plots of Gras a 
Function of ka and klpo for N= 4, Conical Configuration, 1 atm 
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Figure 3-6c,d. Three Dimensional and Contour Plots of Gras a . 
Function of ka and kIno for N= 8, Conical Configuration, 1 atm. 
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Figure 3-6e,f. Three Dimensional and Contuur Plots of Gras a 
Function of ka and klpo for N=15, Conical Configuration, 1 atm. 
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a. CONICAL CONFIGURATION BEAM PATTERN 
To demonstrate the improved directivity of the conical 
configuration over that of the planar system, polar plots can be 
devcloped using the axis system as defined in Figure 3-7, where the 
transducer is displaced is in the positive z direction from the plane of 


the resonators. 





Figure 3-7. Axis System for Beam Pattern Plots, Conical 
Configuration. 


Figure 3-8a shows the beam pattern of the conical configuration in 
both the y,z and the z,6=45° plane. As expected, the symmetry in the 
y.z plane about the z axis is lost, with the sound pressure level (SPL) 
in the positive z direction being 8 dB greater than in the negative z 
direction. If one compares the beam pattern of the conical 
configuration to that of the planar configuration, shown in Figure 3-8b, 
one observes that the SPL in the positive z direction for the conical 
configuration is 1 dB greater than for the planar configuration and 


10dB greater than for one transducer alone. Figure 3-8b clearly shows 








that significant directivity can be achieved with only a minor decrease 
in total power gain, Gr. 
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Figure 3-8a. Beam Pattern for the Conical Configuration 


in the y,z Plane and z,@=45° Plane, Pp>= latm. Scale in 
Decibels (dB). 
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Figure 3-8b. Beam Pattern for the Planar and 
Conical Configuration in the y,z Plane, Po= latm. 


Scale in Decibels (dB). 


3. CONFIGURATION DIMENSIONS 
The results of all calculations performed have been presented in 
terms of the dimensionless radius and spacing, ka and klpo. The actual 
dimensions for the radius a and the separation distance Imax for a given 


frequency are obtained using Equation 3-4. 


ar ont no “Ont (3-4) 


Taking ka = ka, = 0.01379 and kino = x for the optimum values of a 


typical system at an ambient pressure of one atmosphere in Equation 
3-4, a plot of the resulting dimensions a and lno as a function of 


frequency can be made, as shown in Figures 3-9a and 3-9b. 
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Figure 3-9a. Optimum Resonator Diameter (a) at 
Resonance Versus Frequency for a System at Po= latm. 
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Figure 3-9b. Optimum Resonator-Transducer Separation 
(Ino) Versus Frequency for Maximum Power Gain. 
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B. ANALYSIS WITH CONSTANT TRANSDUCER ka, AT 50 
ATMOSPHERES 


Two major assumptions that have been made throughout the 
analysis up this point need to now be modified. The first is that the 
ambient pressure is one atmosphere and the second is that the radius 
of the transducer equals that of the resonators (kao = ka,). In the 
following analysis, an ambient pressure of 50 atmospheres 
(corresponding to a depth of 500 meters) will be taken as a more 
realistic value for an operational system. Since the internal acoustic 
impedance of a bubble is a function of depth, the value for ka, of a 
single bubble at resonance will shift from 0.01379 (Pp = latm) to 
0.09754 (Pp = 50atm). Also, the value of kag will no longer be 
considered equal to that of the resonators, but will be held constant at 
various specific values. 

In the previous sections it was shown that little is gained when the 
number of sympathetic resonators is greater than about four; hence 
the following analysis will be restricted to the case N=4. C-nducting 
similar calculations as before, contour plots and surface plots have 
been generated for the cases N = 4, kao = 0.05, 0.10 and 0.20, for 
both the conical configuration (Figures 3-10a through 3-10f) and the 
planar configuration (Figures 3-1la through 3-11f). As in the previous 
analyses where the ambient pressure was taken to be one atmosphere, 
all cases display a local maximum of Gr of almost 2 at around ka = ka, = 
0.09754 and kl =x. These maximum values are tabulated in Table 3-3. 
Also tabulated, though not plotted, are the cases kag = 0.025 and 


0.150. The first entry in table 3-3 for either configuration was 
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calculated holding kag equal to kay, as in previous calculations, vice 
holding kao constant. A review of Table 3-3 and Figures 3-11 and 12 
demonstrates that the value of kag has little effect on the achievable 


power gain. 


conical 


2.8 
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Table 3-3. Values of Gr for Various Values of kao Ranging from 
.025 to .2 and the Parameters Associated with the Local 
Maximum at kl=nx, ka = ka, = O 0.09754 for both the Conical 
and Planar Configurations at Po= 50atm. 
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Figure 3-10a,b. Three Dimensional and Contour Plots of Gras a 
Function of ka and klno for N= 4, Planar Configuration, 50 atm. 
kag = 0.05. 
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Figure 3-10c,d. Three Dimensional and Contour Plots of Gras a 
kag = 0.10. 


Function of ka and klao for N= 4, Planar Configuration, 50 atm. 








4.0 


3.0 






GR 
2 














WY SS. 
& I, Hh Hi 
SUNN KA VAX Xoo 
RMN IAW AO Xoo 
1 ee LRN 


S 
Ss 


a, 


ho/kor 





> 
: nd 1.5 


f 3 


Figure 3-10e,f. Three Dimensional and Contour Plots of Gras a 
Function of ka and klpo for N= 4, Planar Configuration, 50 atm. 


kag = 0.20. 
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Figure 3-11a,b. Three Dimensional and Contour Plots of Gpas a 
Function of ka and klno for N= 4, Conical Configuration, 50 atm. 


kag = 0.05. 
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Figure 3-11c,d. Three Dimensional and Contour Plots of Gras a 


Function of ka and king for N= 4, Conical Configuration, 50 atm. 
kag = 0.10. 
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Figure 3-1le,f. Three Dimensional and Contour Plots of Gras a 
Function of ka and klno for N=15, Conical Configuration, 50 atm. 


kag = 0.20. 





1. BEAM PATTERNS AT 50 ATMOSPHERES 

Since there is no significant variation in gain with the value of kao, 
beam patterns were only developed for the case kag=ka,, shown in 
Figures 3-12 through 3-14. These plots show the maximum SPL in the 
positive z direction for the conical configuration to be 1.5 dB greater 
than for the planar case and 10.5 dB greater than for the transducer in 


absence of the resonators. 


2. CONFIGURATION DIMENSIONS, FIFTY ATMOSPHERES 
The actual dimensions of the resonator radius a and the 
separation distances lo, for a system operated at an ambient pressure 
of 50 atmospheres were calculated using Equation 3-4 Using the 
optimum value ka = ka, = 0.09754 and klpo = 7. Since the values of Ino 
are not depth dependent, Figure 3-9b still applies. Figure 3-15 


displays the relation between optimum resonator radius and frequency 


at 50 atmospheres. 
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Figure 3-12, Beam Pattern for Planar and Conical Configuration 
in the y,z Plane, Po= 50atm. Scale in Decibels (dB). 
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Figure 3-13. Beam Pattern for Planar Configuration in the y,z 
Plane and z,6=45° Plane, Po= 50atm. Scale in Decibels (dB) 
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Figure 3-13. Beam Pattern for Conical 
Configuration in the y,z Plane,and z,6=45° 
Plane, Po= 50atm. Scale in Decibels (dB). 
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Fig. 3-15. Optimum Resonator Diameter (a) at Resonance 
Verses Frequency for a System at Po= 50atm. 
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IV. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

The problem of multiple scattering of an incident plane wave from 
a system of compact resonant scatterers was solved using both 
Tolstoy's method and network analysis. It was shown that these two 
methods of analysis are equivalent and that both predict the existence 
of "quasiresonance”. Network analysis was then applied to a circular 
array of scatterers, where the incident plane wave was replaced by 
spherical waves radiated from an active transducer, located on the axis 
center. The acoustic advantage of this system was quantified by the 
gain in the radiation resistance seen by the transducer in the presence 
of these scatterers compared to a lone transducer. Two configurations 
were analyzed, (1) the transducer in the plane of the resonators and 
(2) the transducer displaced by one-quarter wavelength from the 
plane. Both configurations showed a gain in radiation resistance of 
almost a factor of two and improved directivity over that of a single 


transducer. 


B. CONCLUSIONS 

The numerical analysis of a transducer surrounded by a circular 
array of sympathetic resonators indicates that a gain in radiation 
resistance of about two is obtainable compared to a lone transducer. 


The physical dimensions of such a system operating at an ambient 
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Pee ae I en 


pressure of 50 atmospheres are displayed in Table 4-1 for several 


System Radius (R) 
(m) 
conical] planer 


frequencies of interest. 


Frequency 
(Hz) 


Resonator Radius (a) 


(m) 


















0.329 64.9 | 74.9 
0.0329 6.49 | 7.49 
0.00329 0.649} 0.749 





Table 4-1. Dimensions of Planer and Conical Configurations 
for Various Frequencies, P, = 50 atm., R is the radius of the 
array. 


Although the calculations presented here were for spherical 
bubbles, the results also apply to more realistic devices (e.g.thin 
shelled cavities, balloons, etc) as long as they are compact (ka<<]). 
M. Strasberg [Ref. 7] showed that the shape of a resonant body has 
little effect on its fundamental resonance frequency. For example, he 
calculated that an oblate spheroid with a ratio of lengths of major to 
minor axes of four, the fundamental resonance frequency differs by 
only a factor of 1.08 from that of a sphere with the same volume. 
Therefore, any deformation of the resonators from a true spherical 
shape, for example due to forces such as buoyancy, should have little 


affect on system performance. 
C. RECOMMENDATIONS 


The results of the research reported in this thesis are very 


encouraging. It is recommended that follow-on work include the 
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design, construction, and testing of an array to verify the theoretical 
predictions presented. 

Theoretical analysis remains to be conducted in at least four areas. 
First, 1n all of the calculations performed, the sympathetic resonators 
were assumed to be gas-filled air bubbles executing radial oscillations 
only. A more realistic model of a sympathetic resonator may be 
required, based upon the properties of a practical device. It is 
recommended this be investigated. Second, Tolstoy has shown that 
the scattering of an incident plane wave from a linear array of resonant 
scatterers exhibits strong directionality. Perhaps a linear array of 
sympathetic resonators offers better performance gain than the 
circular arrays considered in this research. It is recommended that 
such a configuration be investigated. Third, it was shown that in the 
optimum conical configuration the ring of sympathetic resonators lies 
on a cone with a apex angle of approximately 120 degrees, 
independent of frequency. This suggests that a broad-band, constant 
directivity system could be constructed using a series of concentric 
rings of resonators, each tuned to a slightly different frequency and 
placed in its appropriate location along the surface of such a cone, as 
indicated in Figure 4-1. It is recommended that the properties of 
such a system be investigated. Fourth, in the calculations presented, 
the sympathetic resonators were all assumed to be identical, and their 
locations were assumed to be precisely known, as well as that of the 
transducer. It is recommended that the effects of resonator 


nonuniformity and imprecise location be investigated. 
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Resonator r=a’ 





Resonator r=a 


Figure 4-1. Conical Geometric Configuration Allowing an Expanded 
Frequency Response Using Two Concentric Rings of Resonators of 
Radius a and a’, at a Separation Distance from the Transducer of 1 and 
E. 
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APPENDIX 


PROGRAMING INFORMATION 
Inquires into the computer programs used to develop the system 


parameters and associated graphics should be addressed to: 


Professor S.R. Baker 

Physics Department, Code BA 
Naval Postgraduate School 
Monterey, California 93943 
Tel. (408) 646-2732/2729 
Autovon 878-2732/2729 
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